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ABSTRACT 


A  range  rate  expression  is  derived  starting  from  a  general  relativistic 
Doppler  shift  formula.  This  cxpressslon  differs  from  the  one  presently  used 
by  a  term  F(6,^).  The  inclusion  of  this  term,  which  vanishes  for  near-Earth 
ger retries,  allows  range  rate  determination  to  better  than  1  cm/sec. 
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1 .  INTRODUCTION 


A  general  relativistic  fomula  for  a  two-day  Doppler  shift  has  been 
derived  In  a  previous  report.^  Where  the  location  of  the  ground  based 
transmitter  and  receiver  is  the  same»  >nuch  of  the  general  relativistic 
effect  drops  out  and  the  formula  can  be  written  as 


1  - 


-  K 


_ s 

V  ,R  ,  «  . 

1  + 


(1) 


In  Eq.  1,  f^  and  f^  are  the  received  and  transmitted  frequencies,  respectively, 

and  K  is  the  factor  by  which  the  transmitted  frequency  is  multiplied  at  the 

satellite  before  being  retransmitted  back.  The  geocentric  distances  of  the 

satellite  and  transmitter  are  r^  and  r^,  respectively;  v^,  v^  and  v^  are 

the  velocity  vectors  of  the  transmitter,  receiver  and  satellite  respectively;  and 

2 

c  is  the  velocity  of  light  and  m  =  CM  /c  ,  where  G  is  the  gravitational  constant 
and  M  the  mass  of  the  Earth.  R  and  R  are  defined  as 

C*  C  8 


Even  though  the  receiving  and  transmitting  antennas  are  the  same  physical 
object,  we  still  have  r^  ^  r^  and  v^  ^  v^.  This  is  the  result  of  using  a 
stationary  (Inertial)  coordinate  system  with  the  origin  at  the  center  of  the 
Earth  and  measuring  the  coordinates  of  the  antenna  location  at  two  different 


tines,  i.e.,  the  time  of  transmission  and  the  time  of  reception  of  the  signal. 

Since  the  Earth  is  moving,  the  coordinates  of  the  antenna  location  In  the 

inertial  system  are  a  function  of  time.  In  the  limit  of  r  r  ,  for  a  near- 

Earth  situation,  we  also  have  r^  -►  r  and  v^  v  . 

’  t  r  t  r 


II.  INTERPRETATION  OF  THE  DOPPLER  FORMULA 


3 

Eq.  1  Is  relatlvlstlcally  correct  to  orders  of  (v/c)  included  and  m/r 

2 

is  typically  of  the  other  of  (v/c)  .  The  present  range  rate  determination 
method  Is  depedent  on  Che  reading  of  a  cycle  counter.  Truncation  of  one 
count  produces  an  error  of  the  order  of  0.04  ft /sec  ~  1  cm/sec.  For  this 

3 

reason,  in  the  expression  for  f  /f  ,  terms  proportional  to  (v/c)  ,  which 

^  ^  -5 

bring  corrections  of  the  order  of  10  cm/sec  can  be  neglected;  the  terms 
2 

In  (v/c)  which  bring  corrections  of  the  order  of  I  cm/sec  must  be  retained. 

In  accordance  with  the  arguments  presented  in  Reference  1,  the  distinction 

2 

between  r^  and  r^  and  v^  and  v^  in  the  terms  of  order  (v/c)  can  be  dropped 
without  losing  the  desired  accuracy.  Eq.  1  therefore  becomes 


•  /*♦  I  ^ 

where  range  rate, 

R^  can  be  expanded^  In  terms  of  to  yield 

..  0)  At  r 

Rs  -  -  Rt  +  — R - -  ft  R^  ^  (4) 
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/V  A 

where  R.  is  the  unit  vector  normal  to  R  ,  ui  the  angular  frequency  of  the 
ia  c  r 

Earth  and  At  the  time  elapsed  between  the  transmission  of  a  radio  signal  and 
its  reception  after  retransmission  from  the  satellite.  6  is  the  angle  between 
the  velocity  vector  v^  and  the  displacement  R^. 

In  order  to  develop  a  more  suitable  expression  for  the  terms  in  c 
we  note  that  the  velocity  of  the  transmitter  is 


V 

t 


w  X  r. 


(5) 


and  that  the  location  of  the  receiver  can  be  written  in  terms  of  r^  and  v^ 
in  the  following  form 


(0  At 


(6) 


Therefore 


V 

r 


c  X  r. 


+  u»  At 


U)  X  v^ 


) 


(7) 


or 


V  "V  +'.)  Atrv, 
r  r  t  J.,t 


(8) 


where  v.  is  p  unit  vector  normal  to  v  .  Equations  4  and  8  can  now  be  used 
■*  ^  *_] 

to  eliminate  v  and  R  in  the  terms  in  c  .  The  result  is 
r  s 


V  /.t  r. 


V  •  K  +  V  •  R  -  V  •  R  -  V  •  R^  ■  -  2r  -f  — - - sin  0  ( v  -  v  )  •  R 

ssttrsst  R^  \st/  jL,t 
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•f  m'  At  r  v.  •  R  +  0 
t  -l-,t  t 


(9) 
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The  quantity  of  order  o)'’  can  be  neglected  In  Eq.  9,  and  substitution  in 
the  expression  for  gives. 


0)  At  rsin  6 


7^-Kil-2-  +  2(-)  +  -  - - - (  V  -  \  •  R,  ^  +  u)  At  r^v,  ^ 

f^  j  c  \c/  c  L  R^  I  8  t;  -i-.t  t  J-,t  t 

(10) 


The  relationship  between  the  vector  quantities  appearing  in  Eq.  10  is 
shown  in  Fig.  1. 

It  is  obvious  in  Fig.  1  that 


V  ■  cos  6  R  +  sin  OR,,. 

C  C  V 


(11) 


and 


\  cosfl-  e)  -  v^  sin  0 


(12) 


V,  •  R  ■  -  cos 


(I-®)- 


■>  sin  6 


(13) 


%  •  *^a.,t  "  ^S  ^ 


(lA) 


Equation  10  can  now  bn  written  in  the  form 


,2  u  At  r 


i  -  2  f  +  2(^)  +  ®  (^8  ^  “  "'t*^**  ®) 


(1)  At  r^sin  6 


(15) 
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Using  th«  fact  that  u)  ■  u((Dr^)  ■  uiv^  and  raarranglng»  we  hava 


f  r  •  • .  2  111  At  V  T  r  V 

-  2  i  +  2(f)  -  — j-i  (i  e  -  ♦)  •!»  e 


(16) 


or 


f 

=  K 


•.2 


1  -  2  f  +  2(f)  - 


(17) 


which  dafinaa  F(6,^). 

2 

Equation  17  dlffara  fron  the  foniula  which  la  being  presently  used  in 
Che  SCLS  range-rate  processing  by  the  added  ten  F(6»6)«  This  ten  is  an 
affect  of  the  finite  tine  needed  for  a  radio  signal  to  complete  the  trip  to 
the  satellite  and  back.  For  near-Earth  orbits.  At  o  and  F(6,(^)  >  o  as  wall. 
Because  of  the  linear  dependence  on  sin  6,  it  can  be  seen  that  F(6,^)  ■  o  whan  6 
"  o.  This  is  the  case  when  the  satellite  is  on  the  horison.  The  reason  for  this 
ia  that  in  this  Unit  ~  ^(^•*^)  must  also  vanish  when  the  satellite  is  on 

a  circular  synchronous  orbit.  This  is  easily  checked  by  verifying  that  under 
auch  conditions,  the  quantity  f(6,^)  is  zero.  f(6,(^)  is  defined  as 

r  r  V 

f(0,^)  -  1  +  ~  sin  0  -  - — -2-  cos  (18) 

\ 

The  relations  between  the  quantities  involved  in  the  case  of  a  circular 
orbit  for  the  satellite  is  Illustrated  in  Fig.  2  where  it  is  shown  that 


r 

s 


+  cos  ^ 


(19) 
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or 


sin  (6  +  ♦)  +  cos  ())  (20) 

The  constraint  that  the  satellite  be  on  a  geosynchronous  orbit  translates 
into  the  following  equality 


(21) 


trtiere  u  is  now  the  angular  velocity  of  both  a  point  on  the  Earth  and  of  the 
satellite.  To  show  that  under  the  constraints  of  Fqs.  20  and  21,  f(6,4>) 
reduces  to  sero,  Eq.  18  can  be  rewritten  in  the  form 

f(e,4»)  ■  (R^v^  +  r^v^  sin  6  -  r^v^  cos  (22) 


Use  of  Eqs.  20  and  21,  yields 

fj.(6»^)  ■  ”*■  fjV^sin  0  -  r^v^sin  (6  +  4>)  cos  -  R^v^cos^  4>j  /R^v^  (23) 

where  the  subscript  c  Indicates  that  the  constraints  have  been  used.  Carrying 
out  the  algebra  gives 

f  (0,4))  -  V  sin  (J)  [ R  sin  4>  -  r^cos  (0  +  4))1  /R  v  E  0  (2A) 

c  tLt  t 

where  the  last  equality  follows  immediately  with  reference  to  Fig.  2. 
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\ 


Fig.  2.  Relationships  Between  the  Vectors  and  Angles  In 
the  Case  of  a  Circular  Satellite  Orbit 
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III.  APPLICATION  TO  RANGE  RATE  DETERMINATION 


Equation  17,  the  expression  for  the  Doppler  shift,  can  be  used  to 
compute  the  range  rate  of  a  satellite.  In  the  present  system,  an  input 
frequency,  f^,  is  used  where 


(25) 


is  fed  to  a  counter  N^.  In  terms  of  Eq.  17,  the  Input  frequency  becomes 


f .  -  ^  T  +  512  -  (l  - 
i  205  A  c  '  c ' 


+  256  F(e,(j)) 


(26) 


where  K  has  been  set  equal  to  256/205. 

The  count  interval  is  preset  to  last  until  the  N^^  counter  has  accumulated 
■  1,048,574  cycles.  A  second  counter,  N2,  is  slaved  to  start  and  stop  with 
the  counter  and  its  input  is  defined  by 


N 


2 


(27) 


where  the  integral  is  carried  over  the  count  Interval  6t.  6t  is  of  the  order 
of  0.5  sec.  It  should  be  noted  that,  in  general,  because  of  the  presence  of  the 
term  F(6,^)  in  Eq.  26,  the  count  is  no  longer  zero  when  r  ■  0. 


Substitution  of  Eq.  26  Into  Eq.  27,  and  rearrangement  of  the  terms  to 
have  and  appear  on  one  side  and  the  satellite  parameters  on  the  other 
side,  yields 


1_ 

6t 


t  +6t 

I  [('-iV '«•♦>!]<“ 

t 

o 


2048N2 


cN„ 


•f  262A0N 


(28) 


Since  N.  is  a  fixed  number  and  N.  is  a  measured  quantity,  the  right  hand  side 
of  Eq.  28  is  determined.  If  we  replace  r  with  <r>  i.e.,  if  we  substitute  the 
average  value  of  the  range  rate  over  the  count  interval  in  place  of  the 
instantaneous  range  rate,  and  use  the  calculated  satellite  ephemeris  data 
for  the  quantities  appearing  in  F(6,(^),  Eq.  28  will  provide  a  way  to  compute 
<f>.  It  should  be  noted  that  since  the  term  F(6^^)  j  is  of  the  order  of  r^/c  « 
the  knowledge  of  the  quantities  appearing  in  F(6,^)  need  not  be  extremely  accurate. 

In  partluclar,  F(6,(^)  can  be  evaluated  at  any  time  during  the  count  interval  and 

•  • 

can  be  treated  as  a  constant.  While  the  replacement  r  <r>  is  strictly  correct 
t  +6t  t  +6t 

f- 

for  the  term  I  r  dt,  it  is  not  rigorous  for  the  term  I  r  /c  dt  and  in  effect  it 

•4  -4 


o 

amounts  to  the  replacement 


0 


-  dt 
c 


-t  +6t 


<r>  -  dt 
c 


(29) 


This  is  permissible  because  the  quantity  on  the  right  side  of  Eq.  28  is  very 
nearly  the  range  rate. 


-10- 


•  W 


Using  the  definition 


RRN  = 


cN, 


2048N2  +  26240Nj^ 


(30) 


and  the  replacements  suggested  above,  the  following  quadratic  equation  for 
<r>  is  obtained 


<T>^  -  c  <r>  +  c  RRN  -  F(e 


RRN  -  F(e,()))  I ]  -  0  . 


(31) 


The  physically  acceptable  solution  of  Eq.  31  is 


<r> 


RRN  -  F(e,(j))  I 


1  -  [^RRN  -  F(e,.^)  ij  ji 
l3 


(32) 


where  terms  of  the  order  of  jRRN  -  F(e,4>)  jj  jc^  <  lO”^  cm/sec  have  been 
neglected. 

All  the  quantities  on  the  right  hand  side  of  Eq.  32  are  either  measured 
or  calculated  and  an  accurate  determination  of  the  average  range  rate  over  a 
count  interval  is  possible.  Non-relativistic  calculations  are  adequate. 


IV.  CONCLUSIONS 

The  two^ay  Doppler  shift  has  been  seen  to  be  independent  of  general 
relativity  effects  to  a  high  degree  of  approximation.  On  the  other  hand, 
corrections  of  the  order  of  magnitude  of  1  cm/sec  have  to  be  made  for  satellites 
at  geosynchronous  altitude.  The  nature  of  the  correction  is  geometrical,  and 
is  needed  because  the  orientation  of  the  radius  vector  Joining  the  transmitting 
station  and  the  satellite  is  a  function  of  time  and  a  finite  length  of  time  elapses 
between  the  instant  in  which  a  radio  signal  is  transmitted  and  received  back  on 
Earth.  The  range  rate  obtainable  through  Eq.  32  is  averaged  over  a  period  of 
about  0,3  sec  and  is  accurate  to  order  of  10~^  cm/sec. 
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